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ABSTRACT 

We report on continued monitoring of the Anomalous X-ray pulsar (AXP) 
IE 1048.1-5937 using the Rossi X-ray Timing Explorer. We confirm that this 
pulsar has exhibited significant pulsed flux variability. The principal features of 
this variability are two pulsed X-ray flares. Both flares lasted several months 
and had well-resolved few-week-long rises. The long rise times of the flares 
are a phenomenon not previously reported for this class of object. The epochs 
of the flare peaks were MJD 52,218.8 ± 4.5 and 52,444.4 ± 7.0. Both flares 
had shorter rise than fall times. The flares had peak fluxes of 2.21 ± 0.16 and 
3.00 ± 0.13 times the quiescent value. We estimate a total 2-10 keV energy 
release of ~ 2.7 x 10^° ergs and ~ 2.8 x 10^^ ergs for the flares, assuming a 
distance of 5 kpc. We also report large (factor of ~12) changes to the pulsar's 
spin-down rate on time scales of weeks to months, shorter than has been reported 
previously. We find marginal evidence for correlation between the flux and spin- 
down rate variability, with probability of nonrandom correlation 6%. We discuss 
the implications of our findings for AXP models. 

Subject headings: pulsars: general — pulsars: individual (IE 1048.1-5937) — X- 
rays: general 



1. Introduction 

Anomalous X-ray pulsars (AXPs) are an exotic manifestation of young neutron stars. 
AXPs are known for their steady, soft X-ray pulsations in the period range of 6-12 s. The 
detection of X-ray bursts from two AXPs has confirmed the common nature of these objects 
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with that of soft gamma repeaters (?Kaspi et al. 2003), another exotic type of young neutron 
star. Both classes of objects are beheved to be magnetars, i.e., powered by the decay of an 
ultrahigh magnetic field that has a magnitude of 10^^-10^^ G on the stellar surface. For 
recent AXP reviews, see Kaspi & Gavriil (2004) and Kaspi (2004). 

One issue in AXP research has been flux stabihty. Historically, two AXPs have been 
reported to be highly flux variable. Oostcrbrock ct al. (1998) collected all published flux 

measurements for AXP IE 1048.1-5937 and concluded that its total fiux varies by as much as 
a factor of 10 between observations spaced by typically 1-2 yr over ~20 yr. Those data were 
from a diverse set of instruments, including imaging and nonimaging telescopes. Similarly, 
flux variability by a factor of greater than 4 was reported for AXP IE 2259.1+586 by Baykal 
& Swank (1996), using data also from a variety of instruments. 

However, long-term Rossi X-Ray Timing Explorer (RXTE) monitoring of the pulsed flux 
of IE 1048.1-5937 by Kaspi et al. (2001) and of IE 2259.1+586 by Gavriil & Kaspi (2002) 
using a single instrument and set of analysis software showed no evidence to support such 
large variability.^ Also, Tiengo et al. (2002), following a short XMM-Newton observation 
of IE 1048.1-5937, compared the observed flux with those measured by two other imaging 
instruments, ASCA and BeppoSAX. They found that in the three observations, the total 
flux was steady to within ~30%-50%. They argued that the nonimaging detections included 
in the Oostcrbrock ct al. (1998) analysis may have been contaminated by other sources in 
the instruments' fields of view; in particular, the bright and variable X-ray source r] Carina 
lies only 38' away. 

A possible solution to this puzzle came with the discovery of a large (greater than 10 
times) long-lived flux enhancement from IE 2259.1+586 at the time of a major outburst in 
2002 June 18. This event was accompanied by many other radiative changes as well as by a 
large rotational spin- up (Kaspi et al. 2003; Woods ct al. 2004). This suggests that past flux 
variability reported in AXPs could be attributed to similar outbursts that went undetected. 

We report here, using data from our continuing RXTE monitoring program, the dis- 
covery of significant pulsed flux variability in IE 1048.1-5937. This variability is mainly 
characterized by two long-lived pulsed flux flares, having well-resolved rises a few weeks 
long. These are unlike any previously seen flux enhancements in AXPs and SGRs and thus 
likely represent a distinct physical phenomenon. We find no evidence for any major associ- 
ated bursting behavior. We also report large variations in the spin-down torque on timescales 
of a few weeks/months. We find only a marginal correlation between the fiux and torque 
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variations. We argue that this poses another significant challenge to any disk- accretion model 
for AXPs, but is not inconsistent with the magnetar model. 

2. Analysis and Results 

All observations reported here were obtained with the Proportional Counter Array (PCA 
Jahoda et al. 1996) on board RXTE. The timing observations described below are a continu- 
ation of those reported by Kaspi et al. (2001). We refer the reader to that paper for details of 
the analysis procedure. This RXTE monitoring program has shown that, in general, AXPs 
have sufficient stability for phase coherent timing (see Kaspi & Gavriil 2004, for a review). 
IE 1048.1-5937 is an exception. For this pulsar, we have achieved phase-coherent timing 
only over relatively short data spans. In 2002 March, we adopted the strategy of observing 
this source every week with three short (~2 ks) observations. These closely spaced obser- 
vations allow us to measure the spin frequency with high precision weekly without phase 
connecting over long baselines. This therefore allows us to determine the spin-down rate 
with interesting precision on timescales of a few weeks. Figure lA shows the long-term spin 
history of IE 1048.1-5937 as measured by RXTE. 

Figure 2A shows the spin-down rate z> as a function of time over the interval for which 
we can make this measurement. Plotted values of z> were calculated by measuring the slopes 
of each 5 adjacent values of v. Note how v clearly varies greatly during our observations, on 
all timescales to which we arc sensitive. From MJD 52,400 to MJD 52,620 i> had changed by 
a factor of ~12. During the ~120 day interval from MJD 52,620 through MJD 52,740, z> was 
a factor of ~4 larger than the long-term average spin-down ((z>) = —6.48 x 10^^^ Hz s~^). 
This was followed by an abrupt decrease in magnitude by a factor of ~2, which was not 
resolved, and by subsequent additional variations. At no time did we observe any episode of 
spin-up. 

We also monitor the pulsed flux of this source. In this analysis, data from each ob- 
serving epoch were also folded at the optimal pulse period. We calculated the rms pulsed 
flux using the method described by Woods et al. (2004).^ Given IE 1048.1-5937's highly 
sinusoidal pulse profile we used only the first two harmonics to calculate the pulsed flux. 
This method of measuring flux is different from the one used in Kaspi et al. (2001) which 
involved fitting a spectral model to extract a pulsed flux in cgs units. Given the short length 
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of the observations, fitting a spectral model to the individual observations was not practical. 
The pointing for two observations was slightly off-source so we had to correct for reduced 
collimator response. The pointing was on-source for all other observations. Figure IB shows 
our pulsed flux time series in the 2-10 keV band. Pulsed flux time series in the 2-4 and 
4-6 keV bands look similar. 

The pulsed flux time series clearly has significant structure. The most obvious features 
are two long-lived flares. The first fiare was smaller and shorter lived than the second. The 
latter clearly displayed significant structure in its decay In estimating the following fiare 
properties, we define the first fiare as having occurred between MJDs 52,198 and 52,318 and 
the second having started on MJD 52,386, and we take its end to be our last observation on 
MJD 53,030, although it clearly has not yet ended (see Fig. 2). We estimate that the first 
flare had a peak flux of 2.21 ± 0.16 times the quiescent pulse flux, with the peak occurring 
at MJD 52, 218.8 ± 4.5. Its rise time was 20.8 ± 4.5 days, and its faU time 98.9 ± 4.5 days. 
The second flare peak was on MJD 52,444.4±7.0, and had a peak value of 3.00 ± 0.13 times 
the quiescent pulsed flux. Its rise time was 58.3 ± 7.0 days, and its fall time is greater 
than 586 days. We estimate 2-10 keV fiuences of (111 ± 12) x 10^ counts PCU~^ and 
(1136 ± 38) X 10"* counts PCU~^ for the first and second fiare, respectively. Tiengo et al. 
(2002) measured a total fiux in the 2-10 keV energy range of ~ 5 x 10"^^ ergs cm~^ s~^ 
and a pulsed fraction of ~ 94% (for energies >2 keV) from XMM-Newton observations of 
IE 1048.1-5937. This information, along with our measured quiescent pulsed flux, allows us 
to scale our fluences to estimate the total energy released in each flare. Assuming a distance 
of 5 kpc (see discussion in Ozel et al. 2001), we flnd a total energy release of ~ 2.7 x 10^° 
ergs for the flrst flare and ~ 2.8 x 10^^ ergs for the second flare, both in the 2-10 keV band. 

Although we clearly detect both large flux variations and large changes in the spin-down 
rate, the correlation between the two is marginal. The Spearman rank order correlation 
coefficient rs = 0.28, where indicates no correlation and 1 indicates total correlation. The 
probability of obtaining this value of rs or higher by random chance is 6%. Thus, there 
is marginal evidence of some correlation, equivalent to a > 2(t result. From Figure 2, it is 
clear why any correlation is not strong: for example, z> changes very little during the rise 
of the second fiare, in the interval MJD 52,380-52,420. Also, there is no short-term fiux 
change when z> suddenly reaches its maximum absolute value (near MJD 52,620), nor when 
it abruptly changes by a factor of ~2 around MJD 52,740. 

Hardness ratios (HRs) were measured by comparing the pulsed fiux, as measured by the 
method described above, in the 2-4 keV band to that in the 4-6 keV band. Figure IC shows 
our HR measurements. The mean HR is 0.78. There is evidence for spectral variabihty. The 
reduced of the HR time series is 3.6 for 143 degrees of freedom. However, there is no 
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evidence for any correlation of HR with pulsed flux or torque. Our uncertainties, however, 
are quite large; monitoring observations with an imaging instrument would improve this 
situation. 

Intriguingly, the peak of the first flare was coincident with the epochs during which we 
observed two SGR-like X-ray bursts from the direction of this source in 2001 (Gavriil et al. 
2002, indicated by arrows in Fig. 1). However, we found no other SGR-like bursts in any of 
the remaining data. For a detailed description of our burst-searching algorithm see Gavriil 
et al. (2004). We also searched our folded time series for pulse morphology variations using 
the method detailed by Gavriil & Kaspi (2002). We find no evidence for significant pulse 
profile changes at any epoch in our data set. 

3. Discussion 

The long-lived flux enhancements with well-resolved rises that we have observed in 
IE 1048.1-5937 are very different from previously detected X-ray ffux variations in AXPs 
and SGRs, which show very abrupt rises associated with major outbursts (e.g. Kaspi et al. 
2003; Woods et al. 2004). The long-lived flux decay in those sources has been attributed to 
burst afterglow, which is a cooling of the crust following an impulsive heat injection from 
magnetospheric bursts (Lyubarsky et al. 2002). The much more gradual flux rises that we 
have observed in IE 1048.1-5937 comprise a new phenomenon not yet observed in any other 
AXP, despite several years of careful and frequent RXTE monitoring. These flux variations 
may provide a new diagnostic of the physical origin of the persistent nonthermal emission in 
SGRs and AXPs, since they are not contaminated by burst afterglow. Also interesting are 
the large variations in spin-down rate or torque. Torque variations by nearly a factor of 5 
were already reported from RXTE observations (Kaspi et al. 2001), on timescales of years. 
Here we have shown that the torque can change by at least a factor of ~2 more, and on 
much shorter timescales, namely, a few weeks to months. 

In considering the observed pulsed flux and torque variations, whether they are corre- 
lated is an important issue. Our weekly monitoring of the source unfortunately commenced 
only after most of the first fiare decayed. Prior to that, the monthly observations, taken 
in the form of brief snapshots, did not allow anything about the rotational behavior of the 
source to be determined when phase-coherent timing was not possible. This was the case 
during the first fiare. During the second fiare, the spin frequency was, interestingly, most 
stable during the rise and peak of the fiare. Furthermore, the stable spin-down rate was at a 
lower magnitude than the long-term average. Subsequently, ~ 60 days after the fiux began 
to decay, the rate of spin-down began to increase. Given timing observations during only one 
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flare, it is unclear whether these features are coincidences or not. However, there is no strong 
evidence to support otherwise; similar torque variations were seen in the past and were not 
accompanied by any flaring (see Fig. 1). Significant torque variations unaccompanied by 
severe fiux variability have been noted for IE 1048.1-5937 prior to our RXTE monitoring 
(e.g. Paul et al. 2000). Nevertheless, statistically, the probability that they are uncorrelatcd 
is only 4%; studying Figure 2 suggests that if anything, slope transitions are correlated, if 
not the slopes between transitions. Continued RXTE monitoring will help identify any true 
correlations, particularly if the source exhibits more variability. 

Can the magnctar model explain such behavior? The persistent emission in magnctars 
has a spectrum that is well described by a two-component model, consisting of a blackbody 
plus a hard power-law tail. The thermal component is thought to arise from heat result- 
ing from the active decay of a high internal magnetic field (Thompson & Duncan 1996); 
however, thermal X-ray fiux changes are not expected on as short a time scale as we have 
measured in the absence of major bursts. Thompson et al. (2002) put forth a model in which 
the nonthermal component arises from resonant Compton scattering of thermal photons by 
currents in the magnetosphere. In magnetars, these currents are maintained by magnetic 
stresses acting deep inside its highly conducting interior, where it is assumed that the mag- 
netic field lines are highly twisted. These magnctospheric currents in turn twist the external 
dipolar field in the lesser conducting magnetosphere. These magnetic stresses can lead to 
sudden outbursts or more gradual plastic deformations of the rigid crust, thereby twisting 
the footpoints of the external magnetic field and inducing X-ray luminosity changes. The 
persistent non-thermal emission of AXPs is explained in this model as being generated by 
these currents through magnctospheric Comptonization and surface back-heating (Thomp- 
son & Duncan 1996; Thompson et al. 2002). Changes in X-ray luminosity, spectral hardness, 
and torque have a common physical origin in this model and some correlations arc expected. 
Larger twists correspond to harder persistent X-ray spectra, as is observed, at least when 
comparing the harder SCR spectra to those of the softer AXPs. As noted by Kaspi et al. 
(2001), IE 1048.1-5937's hard photon spectral index (F = 2.9) suggests that it is a transition 
object between the AXPs (F ~3-4) and the SGRs (F =2.2-2.4). Hence, if during the fiares 
IE 1048.1-5937's magnetosphere was twisted to the SGR regime, we expect spectral index 
variations of ~ 0.5. Spectral measurements of such precision are not feasible with our short 
RXTE monitoring observations. 

Decoupling between the torque and the luminosity can be accounted for in the magnetar 
model. According to Thompson et al. (2002) the torque is most sensitive to the current 
fiowing on a relatively narrow bundle of field lines that are anchored close to the magnetic 
pole, and so only a broad correlation in spin-down rate and X-ray luminosity is predicted, and 
in fact is observed for the combined population of SGRs and AXPs (Marsden & White 2001; 
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Thompson et al. 2002). However, for a single source, whether an X-ray luminosity change will 
be accompanied by a torque change depends on where in relation to the magnetic pole the 
source of the enhanced X-rays sits. Similarly, large torque variations, as we have observed, 
may occur in the absence of luminosity changes if the former are a result of changes in the 
currents flowing only in the small polar cap region. 

Note that energetically, the total release in these flares is comparable to, although 
somewhat less than, that in the afterglows seen in SGRs and in AXP IE 2259.1+586 (see 
Woods et al. 2004, for a summary). It easily can be accounted for given the inferred magnetic 
energy of the star. 

Although the magnetar model for AXPs has been spectacularly successful in explain- 
ing their most important phenomenology, the anomalous behavior noted for IE 1048.1-5937 
raises the possibility that perhaps it has a physical nature different from other AXPs. It has 
also been suggested that AXPs might be powered by accretion from fossil disks (Chatterjee 
et al. 2000; Alpar 2001). An increase in luminosity Lx can easily be explained in accretion 
models by an increase in the mass accretion rate M, given that Lx oc M. Transient changes 
in M are perhaps not unreasonable to expect in fossil disk models, given the huge varia- 
tions seen in M of conventional accreting sources. However, in an accretion scenario, we 
expect correlations between luminosity and torque. In conventional disk-fed accreting pul- 
sars undergoing spin-up, one expects z> oc L^^. Such a correlation is seen approximately in 
accreting pulsars, with discrepancies possibly attributable to changed beaming or improper 
measurement of bolomctric luminosities, the former due to pulse profile changes, and the 
latter due to finite bandpasscs (Bildsten et al. 1997). As discussed by Kaspi et al. (2001), 
for a source undergoing regular spin-down as in IE 1048.1-5937, the prediction is less clear; 
the form of the correlation depends on the unknown functional form of the torque. For the 
propeller torque prescription of Chatterjee et al. (2000), we find that Lx oc z>^/^, a much 
stronger correlation than in the conventional spin-up sources. For a change in Lx by a 
factor of ~3 as we have seen in the rise of the second flare, we would expect a simultaneous 
change in z> by greater than 50%, clearly ruled out by our data. Conversely, for the abrupt 
change of z> by a factor of ~2 (near MJD 52740), we expect a change in Lx by a factor of 
~5, definitely not seen. This appears to pose a significant challenge to fossil-disk accretion 
models for IE 1048.1-5937. 

Two infrared observations taken on MJD 52,324 (Israel et al. 2002) and MJD 52372 
(Wang & Chakrabarty 2002) have shown that the IR counterpart of this source is variable. 
However, the pulsed X-ray flux at both those epochs was consistent with the quiescent value. 
Furthermore, even though the X-ray flux has not yet returned to its quiescent value, recent 
observations show that the source's proposed IR counterpart is consistent with the fainter of 
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the two previous observations (Durant & van Kerkwijk 2004). This decoupling between the 
IR and the X-ray flux contrasts with what was observed in AXP IE 2259.1+586, whose IR 
flux increased then decayed in concert with the X-ray flux at the time of its 2002 outburst 
(Kaspi ct al. 2003, Tarn et al. 2004, in preparation). This is puzzhng and suggestive of more 
than one mechanism for producing IR emission in AXPs. 
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Fig. 1. — Spin, flux, and spectral history of IE 1048.1-5937. (a) Observed spin frequencies 
vs. time. The points represent individual frequency measurements. The solid lines represent 
the phase-connected intervals as reported by Kaspi ct al. (2001). The dashed line is the 
long-term average spin-down, (h) Pulsed flux time series in the 2-10 keV band. Arrows 
indicate the times at which the bursts reported by Gavriil et al. (2002) occurred, (c) HR 
as a function of time. The HRs displayed were computed for the pulsed flux in the energy 
range (4-6 keV)/(2-4 keV). 



-10- 



Year 



2002 



2003 



-4 

-8 

-12 

-16 

-20 

-24 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



"T 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



2004 



(A) 



is 



Hi'' 



Ml s 



I I I I j I I I I I I I I I : I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I M H -3.0 

(B) 



i 



\ 



I i: I II II I L_L 



0.0 



-1.0 



-2.0 



52200 



52400 



52600 



52800 



53000 



Time (MJD) 



Fig. 2. — (a) z> vs. time over the interval for which our data allow the measurement. The 
horizontal dashed line denotes the long-term average spin-down rate, (z>). The right-hand 
scale is the fractional difference of z> and the long-term average spin-down rate, (h) Zoom-in 
of the pulsed flux time series in the 2-10 keV band. Vertical dotted lines denote the chosen 
start and end ranges for characterizing the two principal flares. 
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